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Abstract 


This  report  constitutes  the  final  report  to  the  Air  Force  Office  of  Scientific  Research 
under  the  Contract  F49620-93-1-0073  which  commenced  on  December  1,  1993,  for 
the  research  project  on  the  Blockage  and  Scattering  of  Lg  Waves. 

The  objective  of  the  project  is  to  investigate  by  means  of  finite-element  modeling 
of  the  problem  of  blockage  and  scattering  of  Lg  wave  propagation  due  to  lateral 
crustal  heterogeneities,  particularly  from  presumed  explosions  at  test  sites  of  Novaya 
Zemlya  across  the  Barents  Sea  to  ARCESS,  and  across  the  Baltic  shield  to  NORSAR, 
and  from  the  Kola  Peninsula  explosion  across  the  Baltic  shield  to  NORSAR.  We 
specifically  address  the  anti-strain  problem  of  SH  waves  related  to  the  Sn  and  Lg 
waves  so  that  there  would  be  no  Pn,  SVn,  SV  type  of  Lg  waves,  and  their  conversions. 
This  research  project  is  intended  to  provide  a  more  quantitative  understanding  of  the 
blockage  and  scattering  of  the  Lg  wave  propagation  due  to  the  presence  of  a  basin. 

Three  nuclear  explosions  were  particularly  finite-element  modeled  concerning  the 
blockage  and  scattering  of  the  Lg  wave  propagation  due  to  the  presence  of  the  basin 
of  the  Barents  Sea  in  comparison  with  the  Lg  wave  propagation  through  the  Baltic 
continental  shield  without  encountering  the  obstacle  of  the  basin  of  the  Barents  Sea. 
These  three  explosions  are:  (1)  the  December  4,  1988  explosion  from  Novaya  Zemlya 
across  the  Barents  Sea  to  ARCESS  with  an  epicentral  distance  of  1,100  km,  (2) 
the  September  4,  1972  explosion  from  Kola  Peninsula  through  the  Baltic  continental 
shield  to  NORSAR  with  an  epicentral  distance  of  1,310  km.,  and  (3)  the  October 
25,  1984  explosion  from  Novaya  Zemlya  across  the  Barents  Sea  through  the  Baltic 
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continental  shield  to  NORSAR  with  an  epicentral  distance  of  2,300  km. 

Geological  structure  models  for  the  great-circle  paths  of  the  Lg  wave  propagation 
of  the  above  three  presumed  nuclear  explosions  were  constructed  and  designated  as 
(1)  Island  Margin  Model,  (2)  Model  I,  (3)  Model  II,  (4)  Model  III,  and  (5)  Model  IV. 
Island  Margin  Model  is  characterized  by  a  sloping  interface  of  the  crust  and  the  upper 
mantle  with  a  dimension  of  450  km  x  60  km.  Model  I  and  Model  II  extend  from  the 
island  margin  of  Novaya  Zemlya  through  the  Barents  Sea  to  ARCESS.  The  island 
margin  portion  of  the  great-circle  path  is  common  to  both  these  two  models.  The 
Barents  Sea  portion  of  these  two  models  is  modeled  as  an  idealized  symmetrically 
shaped  sedimentary  basin,  which  is  placed  in  the  upper  crust  and  located  at  the 
center.  Under  the  basin.  Model  I  is  assumed  to  be  isostatically  uncompensated  with 
a  horizontal  interface  between  the  lower  crust  and  the  upper  mantle.  Model  II  is 
assumed  to  be  isostatically  compensated  so  that  the  interface  between  the  lower 
crust  and  the  upper  mantle  has  an  anti-root  at  the  Moho.  The  Baltic  shield  alone 
is  designated  as  Model  III.  Model  IV  includes  three  segments,  namely,  the  island 
margin,  the  Barents  Sea  and  the  Baltic  shield.  Model  III  is  now  connected  to  Model 
I  or  Model  II  to  form  Model  IV. 

The  nuclear  source  for  simplicity  is  simulated  by  an  impulsive  explosion  source 
of  the  first  derivative  of  Gaussian  type  of  forcing  function.  This  primary  explosive 
source  is  placed  on  the  leeward  side  of  the  island.  Seismic  waves  so  generated  by 
the  impulsive  source  in  Island  Margin  Model  were  then  used  as  the  input  waves  for 
Model  I  and  Model  II  to  provide  the  continuation  of  the  Lg  wave  propagation  from 
the  explosive  source  to  the  receiving  station. 

Three  center-frequencies  of  the  source,  namely,  0.167,  0.334,  and  0.667  Hz  were 
used  to  investigate  the  effect  of  the  center-frequency  of  the  source  on  the  Lg  wave 
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propagation.  In  order  to  study  the  effects  of  basin  width  and  velocity  of  sediments  on 
the  Lg  wave  propagation,  two  basin  widths,  150  km  and  250  km,  were  used,  moreover, 
three  different  sediment  S-wave  velocities,  2.20,  2.70,  and  3.51  km/sec  were  used  for 
the  basin  for  each  basin  width.  An  average  S-wave  velocity  of  the  granitic  and  basaltic 
crust  3.51  km/sec  was  kept  throughout  all  modelings.  The  effect  of  S-wave  attenuation 
in  terms  of  Qs  on  the  Lg  wave  propagation  is  significant.  A  series  of  Qs  ranging  from 
oo,  150,  100,  and  50  have  been  used  in  the  evaluation  of  the  blockage  of  the  Lg  wave 
propagation. 

For  the  center-frequencies  0.167,  0.334  and  0.667  Hz  of  the  source,  the  effect  of 
frequency  content  on  Lg  is  only  minor.  The  arrivals  both  Sn  and  Lg  are  on  time 
except  the  arrivals  of  Sn  are  more  emergent  for  the  center-frequency  of  0.167  Hz. 
The  amplitudes  of  Sn,  Lg  and  scattered  waves  for  these  three  frequencies  are  nearly 
independent  of  frequency.  As  the  transmission  path  of  Sn  is  principally  through  the 
upper  mantle,  the  Sn  arrival  is  virtually  transparent  of  the  presence  of  the  basin. 

The  Lg  waves  are  well  developed  in  the  absence  of  the  Barents  Sea.  For  a  given 
center-frequency  of  the  source,  the  Lg  wave  propagation  is  dependent  on  the  dimen¬ 
sion  of  basin  width.  The  amplitude  of  the  Lg  wave  decays  more  rapidly  for  a  large 
basin  width  than  for  a  small  basin  width.  Moreover,  the  presence  of  a  basin  causes 
time  delay  of  the  Lg  wave  propagation. 

The  velocity  of  sediments  in  the  basin  plays  a  significant  role  in  the  propagation 
and  blockage  of  the  Lg  waves.  The  velocity  contrast  between  the  sediments  in  the 
basin  and  the  surrounding  crust  influences  the  development  and  appearance  of  the 
Lg  waves.  A  higher  value  of  velocity  contrast  causes  a  lengthening  of  the  Lg  wave 
train. 

Among  all  the  effects,  one  of  the  most  important  factors  is  the  effect  of  attenuation. 
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With  reasonable  Qs  for  the  sediments  in  the  basin,  the  blockage  of  the  Lg  wave  is 
very  much  accentuated.  Therefore,  estimation  of  the  quality  factors  not  only  for  the 
sediments  in  the  basin  but  also  for  the  surrounding  crust  becomes  important. 

Comparisons  of  the  finite-element  modeling  results  of  the  December  4,  1988  No- 
vaya  Zemlya  explosion  to  ARCESS  and  of  the  September  4,  1972  Kola  Peninsula 
explosion  to  NORSAR  with  the  real  data  of  the  seismic  traces  FRAO  and  NOlAO(a), 
respectively  show  good  agreement.  The  Sn  waves  are  nearly  the  same  for  both  the 
explosions.  The  wave  train  of  Lg  is  completely  absent  for  the  former  Novaya  Zemlya 
explosion  and  is  perfectly  developed  for  the  latter  Kola  Peninsula  explosion. 

Comparison  of  the  finite-element  simulated  results  of  the  October  25,  1984  No¬ 
vaya  Zemlya  explosion  with  the  NORSAR  trace  NOl  A0(b)  shows  that  the  time  of  the 
arrival  of  Sn  at  the  NORSAR  station  agrees  with  the  finite-element  modeling  calcu¬ 
lation  perfectly.  The  arrival  of  early  Lg  with  small  amplitude,  which  was  interpreted 
by  Baumgardt  as  a  converted  Lg  at  the  continental  margin  of  the  Baltic  shield,  also 
agrees  the  finite-element  simulated  trace.  The  Lg  waves  as  they  are  propagated  across 
the  Barents  Sea  are  suffered  a  great  deal  of  blockage  by  the  low- velocity  sediments 
in  the  basin.  Although  Lg  to  Sn  conversion  and  Lg  scattering  might  also  have  taken 
place  at  the  Baltic  shield  continental  margin,  such  conversion  and  scattering  did  not 
cause  any  major  changes  in  the  already  small  amplitude  Lg  waves. 

On  the  basis  of  finite-element  modeling  and  the  above  three  direct  comparisons 
between  the  finite-element  modeling  results  and  the  observational  data,  the  blockage 
of  Lg  for  the  propagation  paths  from  Novaya  Zemlya  to  ARCESS  and  to  NORSAR 
is  caused  by  the  presence  of  the  low-velocity  sediments  in  the  Barents  Sea. 


Introduction 


The  monitoring  of  underground  nuclear  tests,  unlike  that  of  atmosphere  and  under¬ 
water  nuclear  tests,  can  be  detected  with  relatively  high  degree  of  confidence.  It 
remains  a  critical  part  of  the  global  verification  system. 

Seismology  provides  a  technical  means  for  monitoring  underground  nuclear  testing 
on  the  basis  of  investigating  seismograms  and  knowing  the  characteristic  general 
properties  of  the  propagation  paths  of  teleseismic  body  and  surface  waves  and  of 
regional  waves.  Teleseismic  body  and  surface  waves  are  recorded  to  a  distance  over 
2,000  km.  Regional  waves  include  Pn,  Pg,  Sn,  and  Lg,  which  are  recorded  to  a 
distance  less  than  2,000  km.  Seismologists  are  now  able  not  only  to  calculate  the 
distance  to  the  seismic  event  and  to  deduce  the  type  of  wave  motion,  but  also  to 
make  yield  estimation  and  event  identification  based  on  Lg  amplitudes. 

The  Lg  waves  have  been  in  continuous  attention  for  nearly  one  half  of  a  century. 
The  problems  of  basic  understanding  of  generation,  propagation,  blockage,  and  scat¬ 
tering  remain  highly  challenged.  In  the  sequel,  we  would  selectively  review  a  few 
contributions  to  the  progress  of  the  investigation  of  Lg  waves,  particularly  pertaining 
to  the  blockage  and  scattering  of  Lg  waves. 

The  Lg  wave  was  first  identified  by  Press  and  Ewing  (1952).  Based  on  their  studies 
of  seismograms  recorded  in  North  America,  they  defined  Lg  to  have  periods  of  0.5 
to  6  seconds,  phase  velocities  between  2.0  and  3.5  km/sec,  and  an  inverse  dispersion 
at  distances  greater  than  20*^.  They  noticed  that  Lg  arrivals  are  purely  continental 
paths  and  are  gradually  disappeared  as  the  oceanic  portion  of  the  path  increases  in 
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length  beyond  100  km.  The  arrivals  are  observed  on  all  three  components  but  are 
larger  on  the  horizontals. 

KnopofF  et  al.  (1975)  and  Panza  and  Calcagnile  (1975)  found  that  a  low  velocity 
channel  is  not  necessary  for  the  existence  of  Lg  waves.  Bouchon  (1982)  used  an 
alternative  method  to  evaluate  the  Green’s  functions  for  elastic  layered  media,  and 
computed  complete  synthetic  seismograms  at  short  distances  (150-350  km)  to  model 
an  earthquake  in  France.  He  concluded  that  the  Lg  arrivals  are  composed  of  multiply 
reflected  post-critical  SH  and  SV  rays.  Both  the  lateral  crustal  heterogeneities  and 
the  irregularities  in  the  shape  of  the  two  major  reflectors  (the  free  surface  and  the 
Moho  discontinuity)  should  most  alfect  the  Lg  wave  propagation.  Cara  et  al.  (1981) 
used  an  array  method  devised  by  Cara  and  Minster  (1981)  to  analyze  data  recorded 
in  northwestern  Sierra  Nevada  and  in  southern  California.  They  concluded  that 
representation  of  Lg  as  a  single  multi-mode  wave  train  could  only  explain  observations 
in  the  initial  portions  of  the  wave  train  with  group  velocity  greater  than  3.2  km/sec, 
while  the  later  portions  of  the  Lg  waveform  are  strongly  affected  by  reflections  and 
diffractions  from  any  discontinuities  in  structure.  Kennett  (1986)  used  a  ray  diagram 
approach  to  explain  qualitatively  the  blockage  of  Lg.  Regan  and  Harkrider  (1989) 
used  a  hybrid  of  propagator  matrix  and  finite-element  methods  to  study  the  effects 
of  continental  margin  on  SH  waves. 

The  efficiency  of  Lg  wave  propagation  highly  depends  on  the  structure  and  lithol¬ 
ogy  of  geological  provinces  along  its  propagation  path.  Ruzaikin  et  al.  (1977)  at¬ 
tributed  the  inefficient  propagation  of  Lg  in  the  Tibetan  Plateau  to  the  variation 
of  crustal  thickness,  where  an  unusual  thickness  of  crust  is  present.  Kennett  et  al. 
(1985)  suggested  that  the  Lg  propagation  across  the  Norwegian  Sea  might  be  suffered 
from  a  blockage  of  regional  crustal  thickening  beneath  a  graben  structure  of  a  basin. 
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Baumgardt  (1985,  1990a),  along  the  same  vein,  attributed  to  the  partial  blockage  of 
Lg  propagation  in  Eurasia  by  the  presence  of  the  Ural  Mountains  from  Semipalatinsk 
test  site  in  eastern  Kazakh  to  ARCESS.  Piwinskii  (1981)  suggested  that  the  missing 
granitic  layer  may  cause  the  blockage  of  Lg  waves.  Chan  and  Mitchell  (1985),  and 
Mitchell  and  Hwang  (1987)  suggested  that  the  reduction  of  Lg  amplitude  might  be 
caused  by  the  attenuation  of  the  thick  low-Q  sediments  as  the  Lg  waves  travel  crossing 
the  Barents  Sea. 

Recently,  Baumgardt  (1990a,  1990b,  1991)  used  the  incoherent  beam  method 
(Baumgardt,  1985)  and  ray  tracing  to  investigate  the  Lg  wave  recorded  at  the  regional 
arrays  NORESS  and  ARCESS,  and  Graefenburg  from  presumed  nuclear  explosions 
on  the  northwestern  Russian  platform  and  at  test  sites  of  Novaya  Zemlya.  Baumgardt 
found  that  the  Lg  waves  are  almost  completely  missing  at  ARCESS  and  very  poorly 
recorded  at  NORESS,  while  ’’the  Lg  energy  of  some  kind  appears  to  get  through  to 
the  Graefenburg  array  at  about  the  time  expected  for  on-time  Lg.”  It  is  apparent  that 
variations  of  Lg  amplitudes  are  closely  associated  with  the  problem  of  blockage  and 
scattering  of  Lg  waves  along  its  propagation  path.  Most  recently,  Cao  and  Muirhead 
(1993)  used  the  finite  difference  method  to  simulate  Lg  propagation  in  a  laterally 
varying  continental  path  and  concluded  that  the  existence  of  a  water  column  and 
topographic  variations  of  sea  bottom  and  Moho  could  be  important  factors  for  the 
severe  Lg  attenuation  along  a  continental  path. 

Upon  Nuttli  (1973)  proposed  a  formula  to  relate  amplitude  and  distance  of  Lg 
waves  as  a  high  mode  wave  traveling  with  minimum  group  velocity,  and  in  turn,  to 
define  a  magnitude  scale  (Mb)  for  source  strength  at  short  periods,  attention  has  been 
further  focused  toward  variations  of  Lg  amplitude  in  event  identification  and  yield 
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Nevertheless,  there  is  still  a  great  deal  of  uncertainties  in  basic  understanding  of 
the  generation  and  propagation  of  Lg  waves  in  a  laterally  heterogeneous  crust  that  is 
relevant  to  yield  estimation  and  event  identification  for  nuclear  test  detection.  One  of 
the  problems  of  Lg  propagation  in  the  earth  crust  has  been  the  blockage  and  scattering 
along  its  propagation  path,  and  the  mechanisms  affecting  its  propagation. 

The  1989  announced  intention  of  the  former  Soviets  to  shift  their  underground 
nuclear  testing  activities  from  the  test  site  near  Semipalatinskin,  western  Kazakh  to 
the  Arctic  island  of  Novaya  Zemlya,  made  the  study  of  regional  Lg  propagation  from 
Novaya  Zemlya  to  the  Scandinavian  arrays  relevant. 

This  research  is  motivated  by  the  works  of  a  number  of  investigators,  who  sug¬ 
gested  that  the  blockage  of  Lg  wave  propagation  in  the  continents  primarily  due  to 
lateral  crustal  heterogeneities,  and  particularly  by  a  series  of  papers  by  Baumgardt 
(1985,  1987,  1990a,  1990b,  and  1991),  addressing  the  problem  of  the  blockage  and 
scattering  of  the  Lg  wave  propagation  from  the  Russian  test  sites  to  the  Scandina¬ 
vian  arrays  due  to  the  presence  of  the  Barents  Sea  Basin.  A  more  quantitative  un¬ 
derstanding  of  the  blockage  and  scattering  of  Lg  propagation  and  their  mechanisms 
from  Novaya  Zemlya  to  the  regional  arrays  may  enhance  the  use  of  Lg  amplitude 
discrimination  as  event  identification  and  yield  estimation  elsewhere. 

In  his  F-K  analyses  of  Lg  azimuths,  Baumgardt  (1990a)  showed  that  for  the  first 
arrival  Lg,  the  Lg  propagation  is  nearly  along  a  great-circle  path  from  source  to 
receiver,  although  the  late-arrivals  of  Lg  can  be  multipathed.  In  this  research,  we 
have  focused  our  attention  on  the  study  of  early  Lg  arrivals,  i.e.,  the  great-circle  path 
Lg  arrivals,  and  on  solving  the  problem  of  anti-plane  strain  by  means  of  finite-element 
modeling  of  SH  wave  propagation  in  geologically  realistic  structural  models.  These 
models  include  the  island  margin  margin  model,  i.e.,  Novaya  Zemlya,  the  Barents 


Sea,  a  deep  sedimentary  basin,  and  the  Baltic  continental  shield.  Although  our  initial 
models  are  still  somewhat  simplified,  they  are  essential  to  gain  basic  understanding 
of  the  blockage  and  scattering  of  Lg  propagation  across  the  laterally  heterogeneous 
crust  geological  province  by  province. 


A  Fast  Computation  Algorithm  Of 
The  Finite  Element  Method  For 
The  Anti-Plane  Strain 
Approximations 


For  the  case  of  anti-plane  strain  in  a  Cartesian  system  of  coordinates  x,  y  and  z, 
the  displacement  components  Ux  and  Uy  vanish  everywhere,  and  =  w{x,y,t)  is 
independent  of  the  coordinate  z.  The  non-vanishing  stresses  are 


dw 


Tyz  =  H 


dw 
dy ' 


The  governing  equation  has  the  form  of  wave  equation: 


(1) 


d'^w  d^w  1  d^w 

'd^  ^  ~d^  ~ ^  '■ 


=  III  p, 


(2) 


where  c  is  the  velocity  of  the  shear  wave,  p,  is  the  rigidity,  p  is  the  mass  density,  and 


S  is  the  applied  source  function. 

In  finite  element  algorithm,  the  functional  formulation  corresponding  to  equation 
(2),  and  the  boundary  conditions,  can  be  written  as 


For  wave  attenuation,  we  may  introduce  a  damping  term  in  equation  (3)  so  that 


S  +  p 


d'^w 

'W 
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is  replaced  by 


S  fi 


,dw 


+  P 


d^w 

'W' 


where  n'  is  related  to  /x,  the  shear  modulus,  and  is  to  be  empirically  determined. 

Equation  (3)  must  be  minimized  with  respect  to  the  displacement  w.  (See  any 
finite  element  text  book,  for  example,  Zienkiewicz  (1977)  for  details.) 

The  simplest  two-dimensional  finite  element  in  space  is  defined  by  nodal  points 
i,  j,  k,  and  straight  line  boundaries  as  shown  in  Figure  l-(a).  Each  nodal  point  has 
one  degree  of  freedom.  The  functions  of  the  displacement  field  w{x,  y,  t)  at  any  point 
within  the  triangle  element  is  approximated  by  the  linear  polynomials 


w  =gi+g2x  +  gzV- 


(4) 


The  three  coefficients  gmim  =  1,2,3)  can  be  evaluated  by  solving  one  set  of 
simultaneous  linear  equations.  These  simultaneous  equations  can  be  obtained  by 
equating  the  values  of  the  displacements  at  the  triangle  nodes  i,  j,  and  k.  Then, 
the  coefficients  g,n  can  be  solved  in  terms  of  the  nodal  coordinates.  Once  these 
coefficients  are  determined,  we  can  write  the  displacement  at  an  arbitrary  point  within 
the  triangle  element  as 

{w}  =  [N]  {wY  =  [Ni,Nj,Nk]  {wi,wj,wkf,  (5) 


where  to®  is  a  column  vector  with  three  nodal  displacements  of  the  triangle  element, 
and  Ni,  Nj,  Nk,  or  [N]  are  the  shape  functions,  which  are  expressed  in  terms  of  the 
nodal  coordinates,  and  have  the  form  (Zienkiewicz,  1977): 


Nfn  —  (^m  T  ^mX  d"  CmJ/)/2^,  ITl  i,  J,  k, 


(6) 


Figure  l-(a). 


A  Single  Triangle  Element 


Figure 


(b) .  TWO  Distinct  Ways  of  a  Quadrilateral  Element  and 
Its  Subdivision  into  Two  Triangles 
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where  A  is  the  area  of  triangle  ij  k: 


^  ^  {{^iVk  -  Vj^k)  +  i^kVi  -  VkXi)  +  {xiDj  -  ViXj)]  ,  (7) 

ai  =  Xjyk  -  VjXk, 

hi  =  Vj  -  yfc,  (8) 

Ci  —  Xk  Xj* 

The  other  coefficients  can  be  obtained  by  cyclic  permutation  of  subscripts  in  the  order 
i,  j,  and  k. 

Taking  the  minimization  of  the  functional  (3)  with  respect  to  the  displacement 
field,  we  have 


dx 

d  {tn} 


=  0. 


(9) 


The  final  system  of  eauation  is 


where 


[M]  {ti;}  +  [C]  {th}  +  [K]  {u;}  +  [5]  =  0, 

(10) 

=  (mass) 

(11) 

=  {damping) 

(12) 

Kij^Y^k^j,  (stiffness) 

(13) 

m]-  =  j  NipNjdxdy, 

(14) 

7  [  dx  dy  dy 

(15) 
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dxdy. 


dx  dx  dy  dy  J 
The  total  displacement  field  within  each  element  can  be  obtained  as 


(16) 


w  =  {biWi  +  bjWj  +  bkWk)/2A, 


(17) 


In  the  finite-element  modeling,  we  adopt  a  four-node  quadrilateral  element  as  the 
basic  element  in  the  global  numbering  system.  The  quadrilateral  element  is  comprised 
by  two  constant  strain  triangles  (2CST).  The  two  distinct  ways  of  a  quadrilateral 
element  and  its  subdivision  into  two  triangles  are  shown  as  in  Figure  l-(b).  The 
stiffness  and  mass  of  each  quadrilateral  element  are  obtained  by  averaging  the  results 
of  the  two  types  of  subdivisions  in  order  to  reduce  the  skew  effect. 

In  equation  (10),  the  global  stiffness  matrix  is  a  NxN  square  matrix,  where  N  is 
the  total  number  of  nodal  equations.  In  order  to  prevent  the  stringent  requirement 
of  the  computer  in-core  storage  of  the  global  stiffness  matrix  for  a  large  sized  model 
involving  long  computational  time,  we  adopt  a  nodal-point-oriented  approach  (Teng, 
1981).  The  zero  matrix  elements  in  the  global  stiffness  matrix  are  discarded  to  achieve 
a  fast  computational  algorithm  (Teng,  1989). 

In  adopting  the  nodal-point-oriented  method  and  taking  advantage  of  the  use  of 
uniform  grid  of  finite  elements,  we  reduce  this  large  NxN  square  matrix  for  global 
stiffness  in  2D  for  a  one-degree-of-freedom  problem  to  a  row  vector  with  only  9  mem¬ 
bers  disregarding  no  matter  how  large  the  size  of  the  finite-element  model  may  be  as 
follows: 


{K}  =  {KuK2.K3,...,K,}. 


(18) 


For  time  integration,  we  use  an  explicit  central  difference  scheme: 


{w{t  +  A<)}  =  {ti;  (t)}  +  {w(i)}  Af, 


{th  {t  +  At)}  =  {th  (f)}  -  [I<][M]-^  {u;(f  +  A<)}  A<. 


Geological  Structural  Models 


Gramberg  (1988)  and  Clarke  and  Rachlin  (1990)  provided  fairly  comprehensive  geo¬ 
logical  maps  of  the  Barents  Sea  and  its  vicinities,  including  Novaya  Zemlya,  the  Kola 
Peninsula,  Cheshkaya  Bay,  just  above  the  Arctic  Circle,  as  shown  in  Figure  2-(a). 
Presumably,  on  the  basis  of  the  geological  information  given,  Baumgardt  (1990b)  has 
constructed  a  NW-SE  geological  cross  section  through  A-B  (Figure  2-(b)),  showing 

(1)  the  absence  of  granitic  layer  under  the  southern  Barents  Sea  Basin,  contrary  to 
the  presence  of  a  granitic  layer  under  the  Pechora  Plate  and  the  Central  Barents  Rise, 

(2)  the  thick  accumulation  of  terrigenous  sediments  under  the  Barents  Sea  Basin,  and 
the  thin  terrigenous  sediments  over  the  Plate  and  the  Rise,  (3)  the  variable  thick¬ 
nesses  of  basaltic  layer  present  under  the  Plate,  the  Basin  and  the  Rise,  and  (4)  the 
bending  of  the  Moho  due  to  the  loading  of  the  thick  accumulation  of  sediments  in  the 
Basin,  and  the  irregularities  of  the  depth  of  the  Moho,  ranging  from  about  30  km  to 
37  km  under  the  Basin,  40  km  under  the  Rise,  and  nearly  50  km  under  the  Plate. 

1  Model  for  Novaya  Zemlya  to  ARCESS: 

Similar  geological  cross  sections  along  the  great-circle  path  from  Novaya  Zemlya  to 
ARCESS,  can  be  approximated  constructed.  From  Novaya  Zemlya  to  the  Scandina¬ 
vian  arrays,  NORSAR,  the  great-circle  path  of  Lg  essentially  traverse  major  geological 
provinces,  (1)  the  island  margin  with  variable  crust  thicknesses  (from  Novaya  Zemlya 
to  the  Barents  Sea),  (2)  the  basin  with  variable  width,  depth  and  sediments  (the 
Barents  Sea  Basin),  (3)  the  shield  margin  (from  the  Barents  Sea  to  the  Baltic  shield). 


16 


and  finally  (4)  the  Baltic  shield. 

The  Barents  Sea  Basin  can  be  characterized  by  variations  of  crustal  thickness  with 
sediments  accumulation  as  much  as  15  km,  and  by  missing  of  a  granitic  layer.  This 
structure  is  based  on  low  magnetic  anomalies  caused  by  the  missing  of  granitic  layer 
and  on  high  in  gravity  anomalies  caused  by  the  elevated  Moho.  As  reported,  the 
granitic  layer  in  the  adjacent  province  is  characterized  by  P-wave  velocities  on  the 
order  of  6.0  to  6.5  km/sec,  and  S-wave  velocities  3.46  to  3.75  km/sec  (Clarke  and 
Rachlin,  1990). 

The  sediments  in  the  Barents  Sea  Basin  assume  P-wave  velocities  from  3.9  to 
5.5  km/sec,  and,  correspondingly,  S-wave  velocities  2.2  to  3.17  km/sec  (Baumgardt, 
1990a,  1991).  The  basaltic  layer  assumes  P-wave  velocities  from  6.5  to  6.8  km/sec 
and  S-wave  velocities  from  3.47  to  3.63  km/sec.  The  upper  mantle  is  assumed  to  be 
normal  with  P-wave  velocities  ranging  from  8.1  to  8.2  km/sec,  and  S-wave  velocities 
of  4.6  to  4.7  km/sec. 

The  density  of  the  sediments  are  taken  to  be  2.3  to  2.5  gm/cm^;  that  of  granitic 
layer  2.67;  that  of  basaltic  layer  3.0  in  the  crust;  and  that  of  ultrabasic  rocks  in  the 
upper  mantle  3.27. 

Two  basic  models,  which  are  referred  to  as  ’’Model  I”  and  ’’Model  II,”  as  shown  in 
Figure  3  were  constructed  to  simulate  the  geological  structures  along  the  path  from 
the  former  USSR  nuclear  test  sites  at  Novaya  Zemlya  through  the  Barents  Sea  to 
ARCESS.  The  island  margin  portion  of  the  path  alone,  which  is  common  to  both 
these  two  basic  models,  is  referred  to  as  ’’Island  Margin  Model.” 
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Sedimentary  Layer 
Sections  of  non-granilic  crust 

Isolines  of  cmstal  thickness  (Km) 


Isoiines  of  Granitlc/Metamorphic 
Layer  thickness  (Km) 


(a) 


Central  Barents 
Rise 


Southern  Barents  Basin 


Pechora  Plate 


^^8  Cah^onate/CaJcareous  Basalt  Layer 

Terrigenous  Sediments  ^  Folded  Rocks 

1%^  Crystalline  Rocks 

(b) 


Figure  2.  (a)  Geologic  Maps  of  Barents  Sea  Basin  (After 
Gramberg,  1988) . 

(b)  NW-SE  Cross  Section  (AB)  Across  the  Basin 
(After  Baumgardt,  1990b) . 


direction  of  Lg  propagation 
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(b)  Model  II  simulating  the  geological  structure  along  the  path 
from  the  former  USSR  nuclear  test  sites  Novya  Zemlya  to  the 
Barents  Sea,  which  is  represented  by  a  sedimentary  basin  as  in 
Figure  2- (a);  in  addition  ,  the  Moho  is  assumed  to  have  been 
uplifted  with  an  identical  anti-shape  of  the  sedimentary  basin 


1.1  Island  Margin  Model: 

Island  Margin  model  is  characterized  by  a  sloping  interface  of  the  crust  and  the  upper 
mantle.  The  model  assumes  a  dimension  of  450  km  x  60  km.  A  simulated  source  is 
located  on  the  leeward  side  of  the  island  of  the  model  (Figure  4-(a)). 

The  right-side  artificially  terminated  boundary  is  placed  on  the  other  side  of  of 
the  island  at  AA’;  the  left-side  artificially  terminated  boundary  is  placed  about  100 
km  from  BB’. 

1.2  Model  I  and  Model  II; 

Model  I  and  Model  II  extend  from  the  island  margin  of  Novaya  Zemlya  through  the 
Barents  Sea  to  ARCESS.  As  the  island  margin  portion  of  the  path  is  common  to  both 
these  two  models,  the  Barents  Sea  portion  of  these  two  basic  models  is  modeled  as 
an  idealized  symmetrically  shaped  sedimentary  basin,  which  is  placed  in  the  upper 
crust  and  located  at  the  center. 

Under  the  basin,  Model  I  is  assumed  to  be  isostatically  uncompensated  so  that  the 
interface  between  the  lower  crust  and  the  upper  mantle  is  horizontal,  and  Model  II  is 
assumed  to  be  isostatically  compensated  and  the  interface  is  intruded  by  an  uplifted 
Moho  with  an  identical  anti-shape  of  the  sedimentary  basin. 

The  S-wave  velocities  of  the  granitic/basaltic  (averaged)  layer  and  the  upper  man¬ 
tle  are  taken  to  be  3.51  km/sec  and  4.7  km/sec,  respectively. 

The  widths  of  the  sedimentary  basin  Model  I  and  II  include  150  km  and  250  km. 

2  Kola  Peninsula  PNE  to  NORSAR: 

The  great-circle  path  from  Kola  Peninsula  PNE  to  NORSAR  virtually  lies  in  the 
geological  province  of  the  Baltic  shield.  Geological  structures  in  the  shield  areas 
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Figure  4-(d).  The  Baltic  Shield  Margin  Model. 


for  example,  the  Canadian  Shield  and  the  Baltic  shield  are  fairly  will  known.  The 
shield  crust  basically  consists  of  granitic  and  basaltic  layers  with  relatively  higher  P- 
and  S-wave  velocities  and  densities  as  compared  with  those  in  a  normal  crust.  The 
velocities  for  the  upper  mantle  ultrabasic  rocks  are  assumed  to  be  those  in  normal 
geological  provinces.  A  geological  structure  for  the  great  path  from  Kola  Pennisula 
PNE  to  NORSAR  would  involve  the  Baltic  shield  margin  and  the  shield  itself.  It 
may  be  appropriate  to  model  this  great-circle  path,  designated  as  Model  III  as  shown 
in  Figure  4-(d),  in  which  the  Baltic  shield  margin  has  a  thin  layer  of  lower  S-wave 
velocity  over  a  thinner  crust,  and  the  crust  of  30  km,  then  gradually  increases  to 
about  40  km  toward  NORSAR. 

3  Model  for  Novaya  Zemlya  to  NORSAR: 

The  model  for  the  great-circle  path  for  the  former  USSR  test  sites  from  Novaya 
Zeinlya  to  NORSAR  includes  the  island  margin,  the  Barents  Sea,  the  Baltic  shield 
margin  and  the  Baltic  shield.  Therefore  ,  the  model  is  the  combination  of  the  models 
for  all  these  provinces, namely,  Model  I  and  Model  III,  which  we  refer  to  as  Model  IV 
as  shown  in  Figure  4-(d).  The  physical  parameters  and  the  geometrical  configuration 
of  Model  IV,  some  of  which  are  given  above,  will  be  discussed  in  more  details  in  the 
later  section  on  comparisons  of  the  finite  element  modeling  results  with  observational 
data. 
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Nuclear  Test  Source  Simulation 


The  nuclear  test  source  for  simplicity  may  be  simulated  by  an  impulsive  explosion 
source  of  the  first  derivative  of  Gaussian  type  forcing  functions.  This  primary  source 
is  placed  on  the  leeward  side  of  Island  Margin  Model  of  Novaya  Zemlya. 

In  view  of  the  predominant  frequency  of  Lg  waves  and  the  availability  of  our 
computer  power,  a  compromise  of  the  frequency  content  of  the  source  was  made. 
Three  center-  frequencies  of  the  simulated  nuclear  source  included  0.167  Hz,  0.334 
Hz,  and  0.667  Hz,  corresponding  to  periods  6  sec,  3  sec,  and  1.5  sec,  respectively  to 
investigate  the  effect  of  the  source  frequency  contents  on  the  Lg  wave  propagation. 

Since  our  finite- element  method  is  an  explicit  one,  it  requires  at  least  10-12  grid 
points  per  wavelength  in  each  spatial  dimension  of  discretization  in  order  to  avoid 
grid  dispersion  and  grid  anisotropy,  and  time  step  to  satisfy  the  CRL  conditions  in 
order  to  achieve  accurate  results.  Therefore,  we  consider  the  lowest  S-wave  velocity  of 
2.20  km/sec  and  the  highest  velocity  of  4.70  km/sec,  and,  for  example,  the  maximum 
dimension  of  the  constructed  finite-element  model  for  Novaya  Zemlya  to  ARCESS 
of  1,450  km  X  60  km,  and  one  degree  of  freedom.  In  this  case  of  the  finite-element 
method  requires  a  total  of  solving  5,273  x  219,  2,637  x  110,  and  1,319  x  55  equations 
for  the  cases  of  0.667  Hz,  0.334  hz,  and  0.167  Hz,  respectively.  The  spatial  dimension 
of  the  element  would  be  less  than  1,098  m,  549  m,  and  275  m,  and  the  time  steps 
0.039  sec,  0.078  sec,  and  0.156  sec,  respectively. 

Seismic  waves  so  generated  by  the  impulsive  source  in  Island  Margin  Model  re¬ 
ceived  at  the  cross  section  BB’  are  then  used  as  the  input  waves  at  BB’  for  Model  I 
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and  Model  II  to  provide  the  continuity  of  the  Lg  wave  propagation  from  an  impulsive 
source  on  Island  Margin  Model  to  the  basin  portion  of  the  model  of  Novaya  Zemlya 
to  ARCESS  in  question  (Figure  4-(b)  or  Figure  4-(c).)  The  input  waves  from  the 
cross  section  CC’  provides  the  continuity  of  the  impulsive  source  for  Model  IV,  for 
the  great-circle  path  from  Novaya  Zemlya  across  the  Barents  Sea,  through  the  Baltic 
shield  to  NORSAR. 
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Finite-Element  Modeling  Various 
Effects  On  Lg 

1  Island  Margin  Model: 

First  we  have  investigated  how  the  Lg  wave  is  propagated  through  the  island  margin 
of  Island  Margin  Model,  the  configuration  of  which  is  shown  in  Figure  4-(a).  The 
sloping  interface  at  the  island  margin  is  about  14  In  the  crust,  the  S-wave  velocity 
as  mentioned  previously,  was  taken  to  be  an  average  velocity  of  granitic  and  basaltic 
layers  3.51  km/sec;  and  that  of  the  upper  mantle  4.7  km/sec. 

Figures  5-(a),  5-(b),  5-(c)  are  the  finite-element  synthetic  seismograms  observed 
at  the  epicentral  distance  from  10  to  250  km  for  an  impulsive  source  of  0.167  Hz, 
0.334  Hz,  and  0.667  Hz,  respectively.  Since  we  have  only  solved  the  anti-plane  strain 
problem  of  SH  waves,  as  expected,  we  only  observed  the  arrivals  of  the  Sn  and  Lg 
waves.  If  we  would  have  solved  the  problems  of  plane  strain,  we  would  have  observed  a 
relatively  sharp  Pn  onset,  an  emergent  arrival  corresponding  to  Sn  onset  disregarding 
anisotropy,  converted  waves  and  the  strong  Lg  arrivals. 

The  arrival  of  Sn  for  an  impulsive  source  with  a  center  frequency  of  0.167  Hz  as 
shown  in  Figure  5-(a)  has  a  comparatively  longer  duration  than  that  with  the  center 
frequencies  of  0.334  and  0.667  Hz  in  the  epicentral  distances  from  410  to  650  km 
as  investigated.  The  Lg  wave  begins  with  the  direct  arrival  of  SH,  which  is  often 
included  in  the  energy  envelop  of  the  Lg  waves. 

In  order  to  illustrate  the  crustal  wave  guide  phenomenon  for  the  energy  penetration 
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Figure  5- (a).  The  finite-element  synthetic  seismogram  of 

Lg  waves  by  an  impulsive  line  source  with  center  frequency 
of  0.167  Hz  as  observed  at  distances  of  10  km  to  250  km 
from  the  source  of  the  Island  Margin  Model 
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(c) .  The  finite-element  syntnetic  seismogram  oi 
aves  by  an  impulsive  line  source  with  center  frequency 
,667  Hz  as  observed  at  distances  of  10  km  to  250  km 
the  source  of  the  Island  Margin  Model 


of  the  Sn  and  Lg  waves  and  to  provide  the  Lg  generating  source  for  subsequent 
continuation  to  large  lateral  dimension  models,  we  have  also  calculated  the  finite- 
element  synthetic  seismograms  at  the  cross  section  of  BB’  as  a  function  of  depth  from 
the  crust  through  the  upper  mantle.  The  epicentral  distance  to  the  cross  section  BB’ 
is  400  km.  Figures  6-(a),  6-(b)  and  6-(c)  are  the  seismic  responses  at  the  cross  section 
BB’  with  an  impulsive  source  corresponding  to  the  center  frequencies  0.167,  0.334, 
and  0.667  Hz,  respectively.  As  the  interface  of  the  crust  and  the  upper  mantle  at 
BB’  for  Island  Margin  Model  is  placed  at  a  depth  of  30  km,  all  the  SH-wave  energy 
is  virtually  trapped  in  the  crust  to  make  it  a  perfect  wave  guide  for  the  Lg  wave 
generation  and  propagation.  However,  for  the  center  frequencies  of  0.167  and  0.334 
Hz,  there  are  significant  energy  leakaged  from  the  crust  into  the  upper  mantle.  For 
the  center  frequency  of  the  source  of  0.667  Hz  (Figure  6-(c)),  the  energy  of  Lg  is  all 
confined  to  the  crust.  For  a  lower  center  frequency  0.167  Hz  such  as  shown  in  Figure 
6- (a),  the  energy  of  Lg  in  the  subsurface  is  transmitted  through  the  upper  mantle  as 
deep  as  a  half  of  the  crustal  thickness  to  a  depth  of  45  km;  that  in  Figure  6-(b)  with 
a  center  frequency  of  0.334  Hz,  the  energy  of  Lg  is  transmitted  through  the  upper 
mantle  to  a  depth  of  about  38  km,  or  a  quarter  of  the  crustal  thickness.  Therefore, 
a  crustal  wave  guide  of  Lg  waves  for  a  transitional  region  of  variable  thickness  of  the 
crust  such  as  for  the  island  margin  model  is  of  frequency  dependence. 

In  Figures  6-(a),(b),  and  (c),  there  are  interesting  events,  which  are  generally  not 
observed  on  the  surface.  These  arrivals  in  the  upper  mantle  are  ahead  of  the  arrival 
of  Sn  in  the  crust  as  a  multiply  reflected  S  wave  in  the  margin  and  then  refracted  and 
transmitted  into  the  upper  mantle.  We  refer  these  arrivals  as  SdR,  where  ’cl’  stands  for 
the  reflected  wave  in  the  crust,  and  ’R’  stands  for  refracted  and  transmitted  through 
the  upper  mantle.  The  amplitude  of  SclR  increses  as  a  funtion  of  depth  as  shown  in 
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Fiure  6- (a).  The  seismic  responses  of  Lg  waves  by  an 
impulsive  source  with  center  frequency  of  0 . 167  Hz  as 
observed  along  the  vertical  section  BB'  of  the  Island 
Margin  Model,  which  are  used  to  excite  the  Basin  Model 
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Fiure  6-(c).  The  seismic  responses  of  Lg  waves  by  an 
impulsive  source  with  center  frequency  of  0.667  Hz  as 
observed  along  the  vertical  section  BB'  of  the  Island 
Margin  Model,  which  are  used  to  excite  the  Basin  Mode 


Figures  6-(a),  (b),  and  (c).  The  average  envelops  of  the  Lg  waves  are  broader  for  the 
source  with  lower  center-frequencies  of  the  impulsive  source. 

One  of  the  significant  findings  for  Island  Margin  Model  is  that  shear  wave  re¬ 
verberation  in  the  averaged  granitic/basaltic  layer  is  clearly  shown  to  produce  an 
expected  group  velocity  for  the  Lg  wave  of  approximately  3.5  km/sec,  and  that  for 
the  Sn  wave  of  4.4  km/sec. 

2  Model  I  and  Model  II: 

We  have  studied  the  following  effects  on  the  Lg  wave  propagation  for  Model  I  and 
Model  11. 

1.  Frequency  content  of  a  simulated  explosive  source,  including  its  center  frequen¬ 
cies,  0.167,  0.334,  and  0.667  Hz. 

2.  Basin  width  and  velocity  of  a  sedimentary  basin  in  Model  I  and  Model  II,  which 
connect  to  Island  Margin  Model  at  the  cross-section  junction  of  continuity  BB’.  Two 
basin  widths,  150  km  and  250  km,  and  three  different  sediment  S-wave  velocities, 
viz.,  2.20,  2.70,  and  3.51  km/sec,  were  used  for  the  basin  for  each  width  of  the  basin. 

The  crust  was  assumed  to  have  an  average  S-wave  velocity  of  the  granitic  and 
basaltic  layers,  namely,  3.51  km/sec.  When  the  sediment  velocity  becomes  3.51 
km/sec,  there  would  be  no  presence  of  the  basin.  The  upper  mantle  assumes  a 
S-wave  velocity  of  4.70  km/sec. 

2.1  Without  Sedimentary  Basin: 

We  begin  with  Model  I  and  Model  II  without  the  presence  of  a  sedimentary  basin, 
i.e.,  the  sedimentary  basin  is  replaced  by  the  material  of  the  crust.  In  the  present 
case,  we  simply  used  an  average  S-wave  of  3.51  km/sec  velocity  for  the  entire  crust. 
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Without  the  presence  of  a  sedimentary  basin,  there  would  be  no  uplifted  Moho  in 
Model  II  so  that  geologically  it  makes  sense,  as  most  of  the  crusts  on  the  earth  are 
isostatically  balanced  particularly  in  a  shield  area.  Therefore,  Model  I  is  exactly 
identical  to  Model  II.  In  the  sequel,  we  use  either  Model  I  or  Model  II  that  means 
the  interface  of  the  lower  crust  and  the  upper  mantle  is  horizontal. 

Figures  7-(a),  7-(b),  and  7-(c)  give  the  finite-element  synthetic  seismograms  for  two 
identical  models.  Model  I  and  Model  II,  in  which  the  sedimentary  basin  is  replaced  by 
the  crustal  material,  with  an  impulsive  source  with  the  center- frequencies  of  0.167, 
0.334,  and  0.667  Hz,  respectively.  In  the  present  case.  Model  I  and  Model  II  are 
simply  a  layer  of  crust  with  a  thickness  of  30  km  over  the  upper  mantle  of  30  km. 
The  source  is  located  on  the  leeward  side  of  Island  Margin  Model.  The  continuation 
part  of  Island  Margin  Model  of  Model  I  or  Model  II  was  excited  by  the  waves  arrived 
at  the  cross  section  BB’  as  the  result  of  an  impulsive  source  located  in  Island  Margin 
Model.  The  range  of  these  figures  is  referred  to  with  reference  to  the  cross  section 
BB’  of  Model  I  or  Model  II.  A  distance  of  400  km  must  be  added  to  the  range  in 
order  to  obtain  the  epicentral  distance  in  Model  I  or  Model  II. 

The  essential  features  of  Sn  and  Lg  and  the  scattered  Sn  and  Lg  for  an  impulsive 
source  with  the  center  frequencies  0.167  Hz,  0.334,  and  0.667  Hz  are  precisely  similar. 
The  emergent  arrivals  corresponding  to  the  Sn  onset  are  clearly  observed.  It  follows 
with  distinct  strong  arrivals  of  Lg  waves. 

3  Model  I: 

3.1  Effects  of  Sediment  Velocity: 

Figures  8-(a)  and  8-(b)  are  the  finite-element  synthetic  seismograms  obtained  from 
Model  I  with  a  basin  width  of  150  km  and  two  different  sediment  velocities,  2.7  km/sec 
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Figure  7-(c).  The  finite-element  synthetic  seismogram  for  the 
Basin  Model  of  Figure  4-(b), (basin  width  =  250  km,  basin  depth 
=  15  km),  with  the  basin  velocity  of  3.51  km/sec,  as  observed 
distances  of  20  km  to  400  km  by  a  source  with  a  center  frequent 
0.667  Hz. 
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and  2.2  km/sec.  The  model  was  excited  by  the  input  waves  given  at  the  cross  section 
BB’  as  the  result  of  an  impulsive  source  with  a  center-frequency  of  0.334  Hz  exploded 
in  Island  Margin  Model  as  shown  in  Figure  6-(b)  as  the  input  to  the  Barents  Sea 
portion  of  the  Novaya  Zemlya  to  ARCESS  path. 

In  both  cases,  in  addition  to  the  arrivals  of  Sn  and  Lg  observed  from  20  to  400  km 
from  BB’,  there  are  scattered  events  from  the  two  limbs  of  the  basin  clearly  shown 
at  the  locations  between  the  input  source  and  the  left  limb  of  the  basin  as  marked 
by  the  dashed  lines.  The  reflections  arriving  later  but  stronger  for  the  basin  velocity 
of  2.2  km/sec  case  (Figure  8-(a))  are  due  to  a  higher  contrast  of  velocity  between 
the  crust  and  the  basin  sediments  in  comparison  with  that  of  the  sediment  velocity 
2.7  km/sec  case  (Figure  8-(b)).  For  a  basin  width  of  150  km,  these  scattered  events 
arrive  at  the  distance  about  200  km  from  BB’;  for  that  of  250  km,  they  arrive  at  the 
distance  about  320  km. 

In  the  locations  closer  to  the  input  source  at  BB’,  the  main  features  as  shown 
in  Figure  8-(a)  and  8-(b)  are  the  emergent  arrivals  of  Sn,  which  are  followed  by  the 
distinct  Lg  arrivals.  While  at  the  locations  adjacent  to  or  beyond  the  left  limb  of  the 
basin,  the  Lg  waves  are  developed  into  a  complicated  long  tremor  or  coda,  with  their 
amplitude  diminishing  with  time.  With  the  presence  of  basin,  the  the  arrivals  of  Lg 
waves  are  delayed. 

3.2  Effects  of  Basin  Width: 

Figures  9-(a)  and  9-(b)  are  the  finite-element  synthetic  seismograms  with  the  config¬ 
urations  of  the  model  identical  to  these  for  Figures  8-(a)  and  8-(b),  except  a  basin 
width  of  250  km  was  used.  Ail  the  essential  wave  characteristics  as  revealed  in  Figures 
8-(a)  and  8-(b)  are  preserved.  However,  for  a  basin  width  of  250  km  in  comparison 
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Figure  10- (a).  The  finite-element  synthetic  seismogram  for  the 
Basin  Model  of  Figure  4-(b), (basin  width  =  250  km,  basin  depth 
=  15  km),  with  the  basin  velocity  of  2.20  km/sec,  as  observed  . 
distances  of  20  km  to  400  km  by  a  source  with  a  center  frequent 
0.167  Hz. 


with  that  of  150  km,  there  are  laterally  scattered  waves  from  the  limbs  of  the  basin. 
The  scattered  waves  are  delayed  by  traveling  an  additional  width  of  100  km  in  the 
basin.  These  scattered  waves  contain  both  the  Sn  and  Lg  arrivals  and  have  nearly 
the  same  characteristics  as  those  of  Sn  and  Lg  waves. 

3.3  Effects  of  Frequency  Content: 

With  the  sediment  velocity  of  2.20  km/sec  and  the  basin  width  of  250  km  as  constant, 
Figures  lO-(a),  lO-(b),  and  lO-(c)  are  the  finite-element  synthetic  seismograms  to 
demonstrate  the  effects  of  the  center-frequencies  of  the  source  0.167,  0.334,  and  0.667 
Hz  on  the  Lg  wave  propagation. 

The  essential  features  of  the  arrivals  of  Sn  and  Lg  as  well  as  the  scattered  events 
due  to  the  presence  of  the  basin  are  very  similar  for  three  different  center-frequencies 
of  the  source.  The  arrivals  of  Sn  are  more  emergent  for  the  center-frequency  0.167  Hz 
than  those  for  that  0.334  and  0.667  Hz.  The  amplitudes  of  both  the  arrivals  of  Sn  and 
Lg  as  well  as  the  scattered  events  are  comparable  for  the  center-frequency  of  0.167 
Hz,  and  those  of  the  Sn  and  Lg  are  higher  than  the  scattered  for  the  center-frequency 
of  either  0.334  Hz  or  0.667  Hz. 

The  model  parameters  used  for  generating  Figures  ll-(a),  lT(b)  and  lT(c)  are 
the  same  as  those  for  Figures  lO-(a),  lO-(b),  and  lO-(c)  except  the  sediment  velocity 
is  now  changed  to  2.7  km/sec.  While  the  other  wave  characteristics  of  Sn  and  the 
scattered  events  due  to  the  presence  of  the  basin  remain  similar,  the  amplitudes  of 
Sn,  Lg,  and  the  scattered  events  are  sharply  decreased  due  to  the  differences  of  the 
sediment  velocity. 
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Fiaure  10- (b).  The  finite-element  synthetic  seismogram  for  the 
Basin  Model  of  Figure  4-(b),  (basin  width  =  250  km,  basin  depth 
=  15  km),  with  the  basin  velocity  of  2.20  km/sec,  as  observed 
distances  of  20  km  to  400  km  by  a  source  with  a  center  frequen 
0.334  Hz. 


BB'  (km) 
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ure  11- (a).  The  finite-element  synthetic  seismogram  for  the 
Basin  Model  of  Figure  4- (b), (basin  width  =  250  km,  basin  depth 
=  15  km),  with  the  basin  velocity  of  2.70  km/sec,  as  observed  , 
distances  of  20  km  to  400  km  by  a  source  with  a  center  frequent 
0.167  Hz. 


h-»-<s:lu  <  CO  0)  a  > 


Figure  11- (b) .  The  finite-element  synthetic  seismogram  for  the 
Basin  Model  of  Figure  4-(b)  ,  (basin  width  =  250  km,  basin  dept 
=  15  km),  with  the  basin  velocity  of  2.70  km/ sec,  as  observed 
distances  of  20  km  to  400  km  by  a  source  with  a  center  freque: 
0.334  Hz. 


Figure  11- (c).  The  finite-element  synthetic  seismogram  for  the 
Basin  Model  of  Figure  4-(b),  (basin  width  =  250  km,  basin  depth 
=  15  km),  with  the  basin  velocity  of  2.70  km/sec,  as  observed 
distances  of  20  km  to  400  km  by  a  source  with  a  center  frequent 
0.667  Hz, 


4  Model  II: 

4.1  Effects  of  Basin  Width: 

Figures  12-(a)  and  12-(b)  show  the  finite-element  synthetic  seismograms  for  Model 
II  with  a  center-frequency  of  the  source  0.334  Hz,  a  sediment  velocity  of  2.7  km/sec 
and  two  basin  widths  of  150  km  and  250  km.  Both  the  Sn  and  Lg  waves  are  well 
developed  in  the  distances  from  20  to  400  km  from  the  cross  section  BB’.  However, 
in  the  case  of  the  basin  width  250  km,  the  amplitude  of  the  Lg  decays  more  rapidly 
than  that  in  the  case  of  the  basin  width  150  km.  The  scattered  waves  due  to  the 
limbs  of  the  basin  are  clearly  observed  in  both  the  basin  widths. 

Figures  13-(a)  and  13-(b)  give  the  finite-element  synthetic  seismograms  for  Model 
II  with  a  sediment  velocity  2.2  km/sec;  otherwise  the  model  is  identical  to  that  used 
in  Figures  12-(a)  and  12-(b).  In  Figures  13-(a)  and  13-(b),  the  scattered  waves  due  to 
the  limbs  of  the  basin  arrive  later  but  stronger  because  of  a  higher  velocity  contrast 
between  the  crust  and  the  sediments  in  the  basin  in  comparison  with  a  sediment 
velocity  2.70  km/sec  as  shown  in  Figures  12-(a)  and  12-(b). 

4.2  Effect  of  Basin  on  the  Blockage  of  Lg 

The  Lg  amplitudes  as  observed  on  the  left  limb  of  the  basin  beyond  the  distance  of  300 
km  from  the  cross  section  BB’  for  the  basin  width  250  km  are  drastically  decreased 
in  comparison  with  those  on  the  right-side  of  the  basin  facing  the  source  and  in  the 
basin  proper.  Figure  12-(b)  shows  clearly  that  the  Lg  waves  have  been  blocked  by 
the  presence  of  the  basin. 

In  the  case  of  the  basin  width  150km,  the  blockage  of  Lg  as  shown  in  Figure  12-(a) 
is  not  as  obvious  as  in  the  case  for  the  basis  width  250  km.  The  amplitudes  of  Lg  for 
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with  the  basin  velocity  of  2.70  km/ sec,  as  obi 
distances  of  20  km  to  400  km  by  a  source  with 
freauency  0.334  Hz. 


53 


fi 

o 

o 


Q) 

^  O 
g  C 
X 

W4J 

03 

V  -H 

ca 

iH 

g  (0 

o  ^ 

U  4^ 
<W  C 
<D 
<D  O 
0-H 
C  P4 
flj  0) 
4J 

03  Q) 
•H  43 
Q 

O 

-P 


d) 

43 

-P 

P 

O 

U 


g 

o 

tH 


g 
(0 
p 
CT* 

O  O 


-P 

<W 

^•H 


■P  P 
(d  0) 
4J 

xs  c 

0)  0) 


o 


o 

•H 

-p 
<D 
jC 
4J 
g 
>1 
03  iH 
O 


CU 

W 

I  I 
^  O 
43 
0)  O 
P  S 
:3 

•H  o 


4-1 

O 


4J 

C 

0) 

g 

0) 


4J 

43 

•H 

0) 

43 


> 

P 
0)  (d 
03 

XX  43 
O  -P 
•H 
03  ^ 
rd 

O 
••  O 
O  P 
0)  35 
03  O 
03 

X  <d 

o  >1 

4:3 


CM 


o 

I  -H 

<D 

-P  I 
44 


43 

O  10 

C  H 


4-1 

O 


I 


•H 

44 

<D 

43 


•  iH 
^  0) 
XX  X5 

^  o 

\  X 

CM 

rH  C 
•H 
CD  03 
P  (d 
P  CQ 

•H 

tu 


[| 

43 

44 
a 
0) 

C 

•H 

03 

<d 

XX 


o 
o 
>1*^ 
44 


o 

44 

g 


g 

03  44 
<d  O 


o  ^ 

CM  n 
n 


>1 

O 

g 

0) 

g 

tr 

0) 

p 

44 


<  c/3  0  o  > 


o 

s  ^ 

O 

o  S 

^  00 

+  ^ 
o 

<u  ko 
^  O  CM 
C  a 


MBBSi 

BSSS5BR||^99BBBS| 

SBBb9H9|||99SB999| 

HIH  BI  wJL  I 


o  u  o<  M 

g  W  3  0)  BJ 

(0  I  I  u 

•H  Tf  o  ja  ^ 

0)  43  o  -P 

W  0)  O  -H 

P  S  U1  > 
0  3  <0 

•H  tJ»<t-l  ® 
^  "H  O  ^  O 


4J  4-1  £  u  O 
C  O  t:»'>  w 

>1  "H  g 
W  .H  Q)  «J 
(1)  43 

+J  <0  O  >1 

CO  ^cg  45 

0)  S  g  • 

B  4«<  C'g  B 
0)  43  4^ 

rH  O  If)  4-1 

0)  C  H  o  O 

I  -I-I  o 

<U  CU  II  >1  ■<t 

p  I  -p 

•H  4J  43  -H  O 

C  (0  P  O  -P 

•H  3  O4  O 

4-1  P  0)  iH  B  N 

CJ  -O  0)  >;  53 

0)  > 

43  <0  C  O  ^ 
Eh  C  -H  C  M  n 
(0  W  "H  n 
(0  W  4H  • 

•  r-t  45  3  00 
-^0)  45 

(0  *0  ^  0) 

w  o  B  <u  0)  o 

I  S  5  43  O  C 

n  P  C  0) 

3  3 

p  tr 


0*000  00000000000 

Csj^^OCOOCM^'OOOOCM^'fJOOO 

r-*t-*t-iT-«T-iCSI<M<N01CJ« 


>  43  P 


J— i-hEUJ  <  tO<DO> 


56 


<D 

O 

<W 

g 

OJ 

U 


XI 

'd 


g 
:  o 

H 


II 


in 
cd 
X 

«w 


4J  M 

m  0) 

4J 
TJ  C 
0)  0) 
srH  >  O 


0  0 

04  P 

g  -- 

d  0)  d 

V)  1 

1  CO 

•H  ^ 

0  X  X 

a) 

X  0  4J 

CO  0) 

0  -H 

^  s:  u  > 

0  d 

d 

•H  0) 

-P  -H 

0^0 

0)  pH 

0  M 

X 

+»  0)  :3 

4-) 

43  W  0 

C  O 

cn'^  w 

>i 

•H  g 

W  rH 

0)  Ai  <0 

0) 

43 

4J 

c 

0) 

g 

0) 

rH 

0) 

I 

Q) 

-P 

•H 

c 

-H 

<W 

0) 

X 


X 

o  in 

c  H 
•H 
04  ii 


o 
•  r^ 

t 

%4 

o 


>1 

g 

X 


X 

04 

0) 


o 
o 
>1-^ 
4J 


a  'O 

'd  c 

C  -H 
d  (A 
d 

rH  X 
0) 


o 

-M 

g 


K 


C 
H 
W  Mh 
d  0 


o  r- 

CNl  VO 


d  'O 
w  o 
I  S 

H  C 
•H 
O  V)  03 

U  d 
d  po  Ii 
tp 
•H 


I 

o 

in 


>1 

o 

c 

0) 

d 

d 


►-iEUJ 


<  CO  0)  O  >  A 


57 


Q)  < 

-P 

U 

O 

<4H 


Jd 


C 
•H 

w^J 

•fc-H 


il 


n 

(d 

p 

Cp^fH 

o  o  a 

g  ^  d 
(/}  I  I 
•H  ^  O 
0)  ^ 

(0  <u  o 
u  s 

d 

CP<w 
•H  O 

PL4 

4J 

I  x: 

o  tr* 


o 


,  B 


O 
•H 
4J 
O 

4J 
C 
>1 
03 

Q) 

4J  TJ 
C  O 

0)  s: 

0  42 
H  O  LO 
<D  C  H 
I  -H 
O  04  II 
-P  I 
•H  4->  42 
C  W  4J 
•H  0  a 
m  P 
u 

0) 

42  TJ 
£2 
<d 


4J  P 
(d  0) 
-P 

c 
0)  a) 
>  o 

p 

0)  (d 
01 

43  42 
O  44 

•H 

01  > 
<d 

0) 
^  o 

o  P 

<D  d 

01  o 

03 

B 

o  >1 
r-  43 


CM 

(M 

o 


B 

M 


o 
o 
>1^ 
44 


CD 

T3 


^  0) 
43  TJ 
w  o 
I  S 


c 

•H 

03 

rd 

43 


O 

44 

B  N 
X  32 


C 
•H 
03  M 

fd  o 


o  ^ 
CM  n 
cn 


0) 

P 

d 

•H 

O4 


£2  O 
•H  m 
03  CM 
(d 

01  11 


>1 

o 

c 

0) 

d 

tr 

0) 

p 


h->-is:uj  <  to  CD  o  > 


the  basic  width  of  150  km  after  the  epicentral  distance  of  620  km,  that  is  beyond  the 
left  limb  of  the  basin,  are  only  slightly  different  from  that  of  less  than  620  km.  The 
Lg  amplitudes  in  this  case  remain  nearly  constant  beyond  the  epicentral  distance  of 
620  km. 

The  blockage  of  Lg  propagation  across  a  sedimentary  basin,  therefore,  is  not  as 
simple  as  one  may  expect.  The  phenomenon  of  the  blockage  of  the  Lg  waves  very 
much  depends  primarily  on  the  dimension  of  the  basin,  and  the  contrast  between 
the  sediment  velocity  and  the  surrounding  crust  velocity.  However,  one  of  the  most 
important  factors  of  Lg  blockage  is  due  to  attenuation  by  the  presence  of  low  Q 
sediments  in  the  basin.  We  will  address  the  Lg  attenuation  and  its  blockage  in  the 
sequel. 

4.3  Effects  of  Center- Frequency: 

Figures  14-(a),  14-(b)  and  14-(c)  are  the  finite-element  synthetic  seismograms  for 
Model  II  with  a  basin  width  250  km,  and  a  sediment  velocity  2.7  km/sec  held  constant 
but  varying  the  center-frequency  of  the  source,  0.167,  0.334,  and  0.667  Hz.  Figures 
15-(a),  15-(b)  amd  15-(c)  are  the  same  as  Figures  14-(a),  14-(b),  and  14-(c)  except  the 
sediment  velocity  is  now  2.2  km/sec.  While  a  pronounced  appearance  of  emergence  of 
both  the  Sn  and  Lg  waves  as  well  as  the  scattered  waves  is  observed  for  a  lower  center- 
frequency  of  the  source  0.167  Hz.  As  the  center-frequency  of  the  source  becomes 
higher,  both  the  Sn  and  Lg  as  well  as  the  scattered  waves  have  richer  high  frequency 
contents.  It  is  clear  that  the  higher  the  center-frequency  of  the  source  becomes,  the 
more  the  waves  are  attenuated. 


58 


59 


■■■■■llllMIW 


ft, 


e  -p 

«J  <H  T!  C 
P  Q)  0) 

&>^rH  >  O 
O  U  O.  P 

B  w  3  0)  «} 

U  I  I  U 
•H  ■«i<  O  ,Q  £ 
fl)  X!  O  4-> 

MOO  -H 

P  S  M  > 


uvtmiMiii 


<1-1  4::  M  O 
C  O  D''^  M 
>1  -H  B 
M  -H  O  (0 

O  jG 

4J  >0  O  >1 

CO  A 

OSS* 

g  CM  B 

ox:  M 

r-l  O  in  <1-1 

O  C  H  O  o 


•H  4J  J3  -H  O 
C  W  X>  O  -P 
•H  c  o<  O 
<1-1  P  O  <-l  B  N 
O  -O  O  X  K 
O  > 

XI  TJ  C  Of' 

E-i  C  -H  C  CM  VO 
(0  M  -H  VO 
fl)  M  <H  • 
.  rH  n)  00 


0000000000000000 

c\1'<4'<J3000CM''4’VOOOOCM''TVDCOO 

t-*T-4T-<rHt-<CMC»JC'JCMCMCO 


w  o  B  O  O  O 

I  gAix:  o  c 

■M*  -P  C  O 

I-I  C  o  fl)  c 
•H  in  iS  4->  tJ< 

O  M  CM  -P  M  O 

p  fl)  -H  tH  P 

c  m  II  >  T)  <*-! 


I— 1-4EUJ  <  to  0)  O  > 


60 


0) 

U 

O 

4-1 


s 

o 

H 


II 


fO  4-4  'd 
P  -H  0) 

O  U  Oi  M 
d  0) 
I  I  u 
'I*  O  A 

Xl  o 
0)  o 
^  S  u) 
9  m 
CptM 
r-|  O  ^ 


g 

CO 

•H 

O 

CO 

o 


4J 

0) 

s: 

-p 

c 

>1 

CO 

4J 

g 

Q) 

g 

0) 

O 

I 


(1) 

o 

o 

g 

'H 


4J  0) 

x:  CO 

•H  S 
O  M 
X! 

o 
«-oa 
g  • 

X  rvJ 


to  4^ 
H  O 


P 

0) 

4-) 

g 

0) 

O 

d 

43 

4-> 

•H 

0) 

O 

u 

g 

o 

CO 

d 

g 

4^ 


Q)  04 
4J  I 


•H 

g 

•H 

U 

0) 

Jg 


x: 

-p 


o 
o 
>1^ 
44 


04  O 


•  iH 
^  0) 
d 'd 
w  o 
I  S 
to 

g 

•H 


0 

g 

•p 

CO 

(0 

43 


O 

I 


N 

s 


g 

CO  4-1 
(d  O 


O  C^ 
eg  VO 
H 


0) 

P 

g 

tP 

•H 

04 


I 

o 
to 

CO  04 
d 

n  11 


>1 

o 

g 

(U 

g 

tr 

0) 

p 

4-i 


H- ►-Hsruj  <  CO  0)  o  >  A 


Mlljj 

■mill' 


IHM 


i  MM  MM  Hill  I 

■■■Biii| 


C  O 

Jb4  “H  H 

O  10 

*H  m  II  -P 
A  (Q  0) 

6  _  -H 

«J  <H  03  C 
U  v-H  0)  0) 
>  O 

O  U  a  P 

B  W  3  0)  flj 

(0  I  I  10 
•H  ^  O  A 
0)  J3  O  4J 
to  (U  O  'H 
P  S  (0  > 
O  3  Rj 
•H  O'**-!  0) 

4J  •r^  O  -  O 
(U  b  UP 

jC  4J  Q)  3 
4J  (M  £  CO  O 
C  O  O''^  to 

>1  'H  6 

to  iH  0  44  R3 


>1  'H 
to  iH  0 
0)  4:! 
-P  t3 
C  O  * 
U  S  S 


NHHillll 


00000000000000 

’^VOCOOCsJ^'OCOOW^OCOO 

r-'T-<'i-iT-*r-*CNlCNlC>JCslCvlM 


p  t3  o  >1 
CO  -  OJ  43 

a)  g  s  ♦ 

g  44  M  g 
0)  43  44 

H  0  in  <W 

0)  C  H  O  O 

I  -P  o 

Q)  (X  II 
P  I  4J 
H  +)  43  -H  O 
C  to  -P  U  -P 

r(  s  a  o  ^  • 

IW  P  0)  iH  g  N 

oxj  (a  XX 
<u  > 
c  TJ  (3  o  ^ 

Eh  C  -P  C  CJ  n 

(tj  to  ’H  n 

(d  to  «M  • 

•  tH  43  Id  o  o 

^  0)  43 


I  g  44  43  O 


-P  in  43  4J  O’ 

0)  (0  cs  -P  to  0) 

P  (d  -P  ’P  P 

s  m  II  5  *0  ip 


H-  LU  <  CO  0)  O  > 


62 


ocioooooooooooooo 

CsJ'^^DCX>OCSI'^'XDCX)OCSI'4’UDCOO 
T-iTHT-Hr-4T-HCSICSlC^CMCMCn 


0 

4J 

u 

o 

(M 


s: 

4J  - 
•H  g 

C  o 

H 


II 


e 
(0 
u 

CP 

o  o 


0) 

<0 

W4J 


O 
•H 
4J 

0)  PL4 

x: 

<4-1 

C  0 
>1 

W  H 


a 

w  ;3 

I  I 

^  o 

Q)  o 
CP<w 

o 


-p 

c 

Q) 

g 

0) 

rH 

0) 

I 

0) 

4-) 

•H 

C 

•H 

<4-1 

a) 

43 

EH 


4J  P 
rd  0) 

TJ  C 
0)  0) 
>  O 
P 

0)  (d 

cn 

43  43 
O  44 

•H 

(0 

(d 

0) 
^  o 

o  u 

0)  :3 

U]  o 
w 
g 

44  <d 


43 

0  in 

C  H 


O 

■CNJ 

t 

<4-1 

o 


>« 

43 

g 

44 


0) 

0  TJ 
w  o 
I  S 
in 

H  C 
•H 


0) 

P 

:3 

CP 

•H 


II 

43 

4J 

04 

0) 

T3 

(3 

•H 

in 

CIS 

43 


o 
o 

44 


O 

44 

44 


N 

33 


C 
H 
0}  <4-4 
CIS  O 


O  t-- 
OJ  VO 
VO 


W  CM 

(d 

P3  li 


o 

c 

0) 

p 

CP 

0) 

<w 


I— H-i2:LU  <  C>0  0)  o  > 


4.4  Effects  of  Sediment  Velocity  in  the  Basin 

The  effects  of  the  velocity  of  sedimentary  basin  on  both  the  development  and  appear¬ 
ance  of  Sn,  and  Lg  are  substantive.  When  the  contrast  of  the  velocities  between  the 
basin  and  the  surrounding  granitic/basaltic  crust  increases,  the  wave  train  of  the  Lg 
is  lengthened  drastically. 

4.5  Effects  of  Uplifted  Moho  of  Model  II  on  Lg: 

Although  an  uplifted  Moho  appears  to  have  only  a  minor  effect  on  the  Lg  propagation 
(See,  for  example.  Figures  13-(a)  amd  13-(b)),  an  uplifted  Moho  associated  with  a 
low- velocity  sediment  basin  does  have  some  effects  on  Lg  propagation. 

(1)  The  larger  the  velocity  contrast  between  the  basin  sediments  and  the  surround¬ 
ing  crust  is,  the  stronger  is  scattering,  and  the  longer  are  the  delay  and  attenuation 
of  the  Lg  waves. 

(2)  The  wider  the  basin  width  is,  the  longer  is  the  delay  of  the  Lg  waves. 
Whether  it  is  Model  I  or  Model  II  without  or  with  the  presence  of  an  uplifted  Moho. 
There  are  several  wave-packet-like  arrivals  observed  beyond  the  distance  of  the  left 
limb  of  the  basin.  The  ones  followed  Sn  may  be  interpreted  as  a  conversion  of  the 
Sn-to-Lg  wave  from  the  interface. 

5  Attenuation  and  Blockage  of  the  Lg  wave  Prop¬ 
agation: 

So  far  we  have  not  considered  the  effect  of  sediments  on  attenuation  of  the  Lg  waves. 
Chan  and  Mitchell  (1985),  and  Mitchell  and  Hwang  (1987)  pointed  out  that  the 
reduction  of  Lg  amplitude  might  be  caused  by  the  attenuation  of  the  thick  low-Q 
sediments  as  the  Lg  waves  travel  crossing  the  Barents  Sea. 
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Now  we  introduce  attenuation  in  the  sediments  of  the  basin  with  quality  factors 
Qs  in  finite-element  modeling.  Consider  a  series  of  basin  models  with  different  atten¬ 
uation  values  of  quality  factors,  Qs  =  oo,  150,  100,  50,  with  a  fixed  basin  velocity  of 
2.2  km/sec  and  a  fixed  basin  width  of  250  km. 

Figures  16-(a),  16-(b),  16-(c),  and  16-(d)  are  the  finite-element  synthetic  seismo¬ 
grams  observed  at  the  epicentral  distances  from  620  to  800  km  for  Model  I  with  a 
basin  width  250  km,  sediment  velocity  2.20  km/sec,  a  center- frequency  of  the  im¬ 
pulsive  source  0.334  Hz,  and  Qs  =  oo,  150,  100,  and  50,  respectively.  Figure  16-(a) 
with  an  attenuation  quality  factor  Qs  =  oo  is  identical  to  Figure  9-(b),  i.e.,  without 
attenuation.  Figures  9-(b)  and  16-(a)  have  shown  that  the  blockage  of  Lg  by  the 
presence  of  the  basin.  However,  Figures  16-(b),  16-(c),  and  16-(d)  show  the  blockage 
of  Lg  even  more  clearly  when  the  attenuation  of  the  Lg  waves  is  considered.  The  Lg 
waves  are  noticeably  attenuated  at  the  epicentral  distances  beyond  the  left  limb  of 
the  basin,  i.e.,  greater  than  700  km. 

Figures  17  gives  a  comparison  of  seismic  waves,  including  Sn,  Lg  and  all  other 
waves,  at  the  same  locations  of  observation  for  the  cases  with  different  Qs  values. 
The  four  traces  in  the  figures  correspond  to  the  ceises  of  Qs  =  oo,  Qs  =  150,  100, 
and  50,  respectively.  Comparisons  of  the  seismic  waves  observed  at  the  distances  of 
250  km,  350  km  and  400  km  from  the  source  are  given  in  17-(a),  17-(b),  and  17-(c), 
respectively.  These  figures  show  that  the  Lg  amplitude  is  highly  attenuated  as  Qs 
decreases  in  value,  from  Trace  1  (Qs  =  oo)  to  Trace  4  (Qs  =  50). 

On  the  other  hand,  the  presence  of  a  low-Q  sedimentary  basin  affects  very  little  on 
the  amplitude  of  Sn,  as  the  Sn  waves  are  refracted  mostly  through  the  upper  mantle. 

For  the  present  models  studied,  we  found  that  thick  low-Q  sediments  cause  serious 
attenuation  and  reduction  of  the  Lg  amplitudes,  but  do  not  cause  any  time  delay  of 
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the  Lg  wave  propagation. 
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of  2.2  km/sec  and  with  Qs  -  150 
of  20  km  to  400  km  from  BB' 
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the  basin  models  of  Figure  4-(b),  with  the  basin  velocity 
of  2.2  km/sec  and  with  Qs  =  100,  as  observed  at  distances 
of  20  km  to  400  km  from  BB' 
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Qs  =  CO 


Qs  =  150 


Qs  =  100 


Figure  17- (a), 
locations 
for  cases 


,e  comparison  of  seismic  "f!! 

Observation  at  distance  of  250  km  from  bo 
different  Qs  values. 
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Qs  =  CO 


Qs  =  150 


Qs  =  100 


Figure  17- (b) .  The  comparison  of  seismic  waves  at  the  same 
locations  of  observation  at  distance  of  350  km  from  BB' 
for  cases  with  different  Qs  values. 
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Qs  =  CO 


Qs  =  150 


Qs  =  100 


Figure  17- (c) . 
locations 
for  cases 


The  comparison  of  seismic  waves  at  the  same 
of  observation  at  distance  of  400  km  from  BB 
with  different  Qs  values. 


Comparisons  Of  Finite-Element 
Modeling  Results  With  Real 
(Observational)  Data 


We  have  learned  through  finite- element  modeling  as  much  about  the  effects  of  island 
margin,  basin,  and  shield  in  terms  of  basin  width,  center-frequency  of  an  impulsive 
source,  and  wave  velocity  of  sediments  in  a  basin,  as  well  as  attenuation  in  terms  of 
quality  factor  on  the  Lg  wave  propagation. 

An  ultimate  test  of  the  goodness  of  finite-element  modeling  is  to  compare  finite- 
element  synthetic  seismograms  based  on  geologically  realistic  structure  models  with 
real  data. 

Among  the  available  data,  it  seems  that  three  Lg  propagation  paths  used  by 
Baumgardt  (1991)  in  his  study  of  the  blockage  of  the  Lg  waves  most  appropriate  to 
make  comparisons  of  finite-element  modeling  with  observational  data. 

(1)  The  ARCESS  FRAO  trace  of  the  December  4,  1988  Novaya  Zemlya  explosion, 
representing  an  oceanic  path.  The  epicentral  distance  from  the  former  USSR  test 
sites  of  Novaya  Zemlya  to  ARCESS  is  about  1 ,100  km.  The  December  4,  1988  Novaya 
Zemlya  explosion  had  a  magnitude  of  5.7. 

(2)  The  NORSAR  NOl  AO(a)  trace  of  the  September  4, 1972  Kola  Peninsula  explosion 
(PNE),  representing  nearly  a  continental  path.  The  epicentral  distance  from  Kola 
Peninsula  to  NORSAR  is  about  1,310  km  slightly  greater  than  that  for  the  above 
Novaya  Zemlya  explosion  to  ARCESS.  The  Kola  Peninsula  explosion  had  a  magnitude 
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of  4.6. 


(3)  The  NORSAR  NOlAO(b)  trace  of  the  October  25,  1984  Novaya  Zemlya  explosion, 
representing  a  mixed  path  from  the  island  margin,  through  the  Barents  Sea  across 
the  Baltic  shield  to  NORSAR.  The  epicentral  distance  of  this  event  of  a  magnitude 
of  5.9  is  2,300  km. 

The  great-circle  paths  of  the  above  three  events  are  shown  in  Figure  18.  The 
playback  traces  of  these  three  events  are  displayed  in  Figure  19,  showing  the  arrivals 
of  Pn,  Sn,  and  Lg.  These  traces  as  pointed  out  by  Baumgardt  (1991)  have  been 
filtered  in  the  0.6  to  3.0  Hz  frequency  band.  Although  the  traces  of  these  events 
have  been  recorded  with  different  short-period  instrument  responses,  these  differences 
should  not  affect  the  relative  amplitude  information  in  a  major  way.  Because  of  the 
difference  in  magnitude  and  instrument  response,  Baumgardt  scaled  the  amplitudes 
of  the  second  trace  of  NORSAR  by  a  factor  of  three  to  the  first  trace  of  ARCESS  so 
that  the  Pn  amplitudes  for  both  traces  are  equalized.  The  third  trace  of  the  Novaya 
Zemlya  to  NORSAR  event  of  October  25,  1984  is  as  playbacked. 

So  far  as  Sn  and  Lg  are  concerned,  comparison  of  the  NORSAR  trace  for  the  Kola 
Peninsula  explosion  with  the  trace  of  ARCESS  for  the  Novaya  Zemlya  event  shows 
noticeably  that  the  wave  train  of  the  Lg  waves  is  well  developed  for  the  ARCESS  trace, 
while  it  is  virtually  missing  for  the  NORSAR  trace.  The  Sn  waves  are  nearly  the  same 
for  both  the  traces.  The  great-circle  path  for  Sn  and  Lg  traveling  from  Novaya  Zemlya 
to  ARCESS  could  be  considered  as  Model  I  or  Model  II  in  the  previous  finite-element 
modeling,  while  that  from  Kola  Peninsula  to  NORSAR  could  be  considered  as  Model 
III,  a  continental  paths  from  the  Baltic  shield  margin  through  the  Baltic  shield. 
Again,  we  have  used  the  method  of  segment  by  segment  finite-element  calculation. 
The  island  margin  portion  of  the  model  is  identical  to  that  used  before,  as  shown  in 
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Figure  18.  Map  showing  comparisons  of  the  propagation 
path  from  the  Kola  Pennisula  to  NORSAR,  the  Novaya 
Zemlya  to  ARCESS  and  from  Novaya  Zemlya  to  NORSAR. 

The  two  paths  from  Novaya  Zemlya  to  ARESS  and  from 
the  Kola  PNE  to  NORSAR  are  nearly  parallel  great-circle 
paths  and  have  almost  identical  epicentral  distance. 


Figure  4-(a),  except  that  now  the  epicentral  distance  to  the  cross  section  BB’  is  500 
km.  The  explosion  source  in  Island  Margin  Model  is  an  impulsive,  transient  source 
of  the  fourth  derivative  of  Gaussian  type  forcing  functions  with  a  center-frequency  of 
0.334  Hz. 

Finite-element  Modeling  of  the  December  4,  1988 
Novaya  Zemlya  Explosion  and  Compared  with  the 
ARCESS  Data: 

On  the  basis  of  Model  I,  synthetic  seismograms  for  the  Novaya  Zemlya  December  4, 
1988  explosion  are  generated  by  means  of  finite-element  modeling  for  an  epicentral 
distance  from  800  to  1,300  km  from  the  explosion  as  shown  in  Figure  4-(a).  The 
Barnets  Sea  is  modeled  as  a  sedimentary  basin.  The  geophysical  parameters  of  the 
model  are  that: 

1.  the  configuration  of  the  basin  is  shown  in  Figure  4-(b),  with  a  basin  width  of 
700  km.;  the  density  of  the  sediments  was  taken  to  be  2.45  gm/cm3;  the  shear-wave 
velocity  was  taken  to  be  2.3  km/sec,  and  Qs  was  100. 

2.  for  the  granitic/basaltic  crust,  the  shear-wave  velocity  was  assumed  to  be  3.58 
km/sec  at  the  top  and  gradually  increased  to  3.75  km/sec  at  the  bottom  of  the  crust; 
its  density  ranged  from  2.67  gm/cm3  to  2.90  gm/cm3,  and  an  averaged  Qs  200. 

3.  for  the  upper  mantle,  the  S  velocity  was  taken  to  be  4.70  km/sec  and  the 
density  3.3  gm/cm3. 

Figure  20-(a)  gives  the  finite-element  synthetic  seismograms  for  the  December  4, 
1988  Novaya  Zemlya  explosion  at  an  epicentral  distance  of  800  to  1,300  km.  The 
arrivals  of  Sn  is  distinct  and  the  signal  of  Sn  attenuates  with  distance.  The  wave 
train  of  Lg,  although  it  also  suffers  attenuation,  disappeared  beyond  the  epicentral 
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Figure  20- (a).  The  finite-element  synthetic  seismogram  of 

Lg  waves  of  the  simulated  model  for  the  ARCESS  trace  of 


H-i-iELJ  <  CO  OJ  O  >  A 


Figure  20- (b).  Comparison  of  the  numerical  result  and  the 
real  data  for  the  ARCESS  FRAO  trace  of  the  December  4 
1988  explosion. 


distance  of  about  1,100  km  on  the  other  side  of  the  Barents  Sea. 

Figure  20-(b)  plots  both  the  ARCESS  trace  and  the  finite-element  modeled  trace 
for  comparison,  showing  the  onset  of  Sn  and  the  absence  of  the  Lg  wave  altogether. 
Agreement  between  the  finite-element  modeling  results  and  the  real  data  is  very 
encouraging. 

Finite-element  Modeling  of  the  September  4,  1972 
Kola  Peninsula  Explosion  and  Compared  with  the 
NORSAR  Data: 

On  the  basis  of  Model  III,  the  Baltic  shield  model,  a  series  of  seismograms  have  been 
generated  at  an  epicentral  distance  of  800  to  1,300  km. 

The  structural  model  of  the  Baltic  shield  as  shown  in  Figure  4-(d)  includes  a 
two-km  thick  low  velocity  layer,  which  has  a  shear-wave  velocity  of  2.6  km/sec  and 
density  of  2.3  gm/cm3,  and  a  Qs  value  of  100.  The  geophysical  parameters  of  the 
granitic/  basaltic  crust  and  the  upper  mantle  are  the  same  as  those  used  in  modeling 
the  above  Novaya  Zemlya  event. 

Figure  21- (a)  shows  the  finite-element  synthetic  seismograms  for  the  Kola  Penin¬ 
sula  event.  The  distinct  features  of  these  seismograms  are  the  prominent  Lg  wave 
onset  across  all  the  epicentral  distances.  There  is  no  question  obout  it  that  the  Lg 
waves  are  propagated  without  any  blockage.  The  Sn  wave  propagation  for  the  Kola 
Peninsula  explosion  has  nearly  the  same  characteristics  of  those  for  the  Novaya  Zemlya 
explosion,  as  they  should  be  because  the  propagation  paths  in  both  the  explosions 
are  refracted  through  the  upper  mantle. 

Figure  21-(b)  shows  the  comparison  of  the  finite-element  trace  with  the  NORSAR 
NOlAO(a)  trace,  giving  the  arrivals  of  the  Sn  and  Lg  clearly.  As  the  center-frequency 
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Figure  21-(a).  The  finite-element  synthetic  seismogram  of 

Lg  waves  for  the  structural  model  of  the  simulated  model 
for  the  NORSAR  NOIAO  of  the  September  4,  1972  Kola 
Peninsula  explosion. 


Figure  21- (b) .  Comparison  of  the  numerical  result  and  the 
real  data  for  the  NORSAR  NOIAO  of  the  September  4. 
1972  Kola  Peninsula  explosion. 


of  the  impulsive  source  in  finite-element  modeling  is  not  high  enough  in  comparison 
with  the  trace  of  the  real  data,  the  frequency  contents  of  both  the  Sn  and  Lg  in 
finite-element  modeling  results  appear  to  be  lower.  However,  all  the  characteristics  of 
the  Sn  and  Lg  waves  through  finite-element  modeling  are  very  similar  to  the  observed 
trace  NOlAO(a). 

Finite- element  modeling  of  the  October  25,  1984 
Novaya  Zemlya  Explosion  and  Compared  with  the 
NORSAR  Data: 

The  above  two  comparisons  of  the  finite-element  modeling  results  with  the  observed 
data  are  for  the  Lg  propagation  for  a  purely  continental  path  and  for  an  oceanic  path. 

From  the  above  two  explosions,  we  have  gained  considerable  physical  insight  to 
the  problem  of  the  blockage  and  scattering  of  the  Lg  wave  as  the  wave  is  propagated 
through  the  island  margin,  and  then  its  continuation  of  across  the  Barents  Sea  through 
the  segment  by  segment  finite-element  calculation.  We  have  assessed  the  influence  of 
the  Barents  Sea  on  the  Lg  waves  in  term  of  frequency,  basin  width,  sediment  velocity, 
and  the  geometry  of  the  crust  and  the  upper  mantle  contact. 

Our  interest  in  this  section  is  naturally  to  investigate  the  influence  of  the  island 
margin,  the  Barents  Sea,  and  the  Baltic  shield  altogether  on  the  Lg  wave  propagation. 
For  this  purpose  we  have  constructed  a  complete  model,  named  Model  IV.  Model  IV 
thus  includes  three  segments,  namely,  the  island  margin,  the  Barents  Sea,  the  Baltic 
shield  margin,  and  the  Baltic  Shield  as  shown  in  Figure  4-(d).  Baltic  Shield  Margin 
Model  is  now  connected  to  Model  I  at  the  cross  section  of  CC’  to  form  Model  IV. 
Adopting  Model  II  with  an  uplifted  Moho  in  Model  IV  would  not  make  any  significant 
changes  on  the  arrival  time  or  amplitude  of  the  Lg  waves,  but  it  would  only  have  a 
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very  minor  difference  in  amplitude  and  arrival  time  of  the  Sn  onset.  Therefore,  in  the 
present  analysis,  Model  IV  adopts  Model  I  with  a  horizontal  Moho. 

All  the  physical  parameters  of  Model  IV  are  exactly  identical  to  Model  I  used  for 
the  comparison  of  the  finite-element  modeling  results  with  the  ARCESS  trace  of  the 
December  4,  1988  Novaya  Zemlya  explosion  except  the  dimensions  of  three  segments 
are  changed  to  that: 

1.  The  termination  cross  section  BB’  of  the  island  margin  portion  of  Model  IV  is 
placed  at  an  epicentral  distance  of  400  km. 

2.  The  termination  cross  section  CC’  of  the  Barents  Sea  portion  of  Model  IV  is 
placed  at  an  epicentral  distance  of  1 ,300  km. 

3.  The  continuation  portion  of  the  cross  section  CC’,  i.e..  Model  III  of  the  Baltic 
shield,  which  connects  to  the  above  two  segments  to  complete  Model  IV. 

Figures  22-(a)  and  22-(b)  show  the  finite-element  synthetic  seismograms  for  the 
continuation  portion  of  Model  IV  which  clearly  show  that  the  Sn  onsets  are  prominent, 
and  the  energies  for  both  early  and  Lg  are  very  weak.  Apparently  the  Lg  waves  have 
already  been  blocked  as  they  are  propagated  across  the  Barents  Sea  Basin. 

Figure  22- (c)  shows  the  comparison  of  the  NORSAR  NOl  A0(b)  of  the  October  25, 
1984  Novaya  Zemlya  explosion  with  the  finite-element  simulated  trace  at  the  same 
epicentral  distance  2,300  km.  The  time  of  the  arrival  of  Sn  at  the  NORSAR  station 
agrees  with  the  finite-element  modeling  calculation  perfectly.  The  arrival  of  early  Lg 
with  small  amplitude,  which  was  interpreted  by  Baumgardt  as  a  converted  Lg  at  the 
continental  margin  of  the  Baltic  shield,  also  agrees  the  finite-element  simulated  trace. 
The  arrival  of  Lg,  which  was  predicted,  had  nearly  the  same  amplitude  as  that  of  Sn. 
Apparently,  the  Lg  waves  as  propagated  across  the  Barents  Sea  have  been  blocked  by 
the  low-velocity  sediments  in  the  basin.  Although  a  Lg  to  Sn  conversion  might  also 


84 


have  taken  place  at  the  continental  margin  of  the  Baltic  shield. 

The  above  three  direct  comparisons  between  the  finite-element  modehng  results 
and  the  real  observational  data  substantiate  that  the  blockage  of  Lg  for  the  propaga¬ 
tion  paths  from  Novaya  Zemlya  to  ARCESS  and  to  NORSAR  is  certainly  caused  by 
the  presence  of  the  low-velocity  sediments  in  the  Barents  Sea. 
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ire  22- (a).  The  finite-element  synthetic  seismogram  of  Lg 
waves  of  the  simulated  model  for  the  October  25,  1984 
Novaya  Zemlya  explosion,  recorded  at  NORSAR  for  epi  central 
distances  1,700  to  2,400  km. 
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Figure  22- (c).  Comparison  of  the  numerical  result  and  the 
real  data  for  the  NORSAR  NOIAO  of  the  October  25, 
1984  Novaya  Zemlya  explosion. 


Conclusions 


1.  Through  extended  finite-element  modeling,  it  is  now  proven  that  large- dimension 
geologically  realistic  anti-plane  strain  problems  of  SH  waves  can  be  efficiently  and 
accurately  modeled  by  means  of  the  method  of  segment  by  segment  calculation. 

2.  Finite-element  modeling  of  the  Lg  wave  propagation  along  the  great-circle 
paths  from  Novaya  Zemlya  to  ARCESS  designated  as  Model  I  and  Model  II,  and 
from  Kola  Pennisula  to  NARSAR  designated  as  Model  III  indicates  that 

a.  For  the  center- frequencies  0.167,  0.334  and  0.667  Hz  of  the  source,  the  effect 
of  frequency  content  on  Lg  is  only  minor.  The  arrivals  of  both  Sn  and  Lg  are  on 
time  except  the  arrivals  of  Sn  are  more  emergent  for  the  center-frequency  of  0.167 
Hz.  The  amplitudes  of  Sn,  Lg  and  scattered  waves  for  those  three  frequencies  are 
nearly  independent  of  frequency.  The  Lg  waves  are  well  developed  in  the  absence  of 
the  large- dimension  of  sedimentary  basin  of  the  Barents  Sea. 

b.  For  a  given  center  frequency  of  the  source,  the  Lg  wave  propagation  is  depen¬ 
dent  on  the  dimension  of  basin  width.  The  amplitude  of  the  Lg  wave  decays  more 
rapidly  for  a  large  basin  width  than  for  a  small  basin  width.  The  presence  of  a  basin 
causes  the  time  delay  of  the  Lg  wave.  As  the  transmission  path  of  Sn  is  principally 
through  the  upper  mantle,  the  Sn  arrival  is  virtually  transparent  of  the  presence  of 
the  basin. 

c.  The  velocity  of  sediments  in  the  basin  plays  a  significant  role  in  the  propagation 
and  blockage  of  the  Lg  waves.  The  velocity  contrast  between  the  sediments  in  the 
basin  and  the  surrounding  crust  influences  the  development  and  appearance  of  the 
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Lg  waves.  A  higher  value  of  this  velocity  contrast  lengthens  the  Lg  wavetrain. 

d.  Among  all  the  effects,  one  of  the  most  important  factors  is  the  effect  of  wave 
attenuation.  With  reasonable  Qs  for  the  sediments  in  the  basin,  the  blockage  of  the  Lg 
wave  by  the  presence  of  the  Barents  Sea  Basin  is  very  much  accentuated.  Therefore, 
estimation  of  the  quality  factors  not  only  for  the  sediments  in  the  basin  but  also  for 
the  surrounding  crust  is  of  considerable  importance. 

3.  Comparisons  of  the  finite-element  modeling  results  of  the  former  USSR  test 
sites  of  Novaya  Zemlya  to  ARCESS  and  to  NORSAR,  and  of  the  Kola  Peninsula  to 
NORSAR  with  the  real  data  of  the  seismic  traces  NOl  A0(a,b)  and  FRAO,  respectively 
show  that 

a.  The  characteristics  of  an  emergence  of  Sn  and  the  expected  time  of  the  arrival 
of  Sn  are  virtually  unchanged  for  all  the  above  events  investigated. 

b.  The  Lg  waves  are  very  well  developed  for  the  September  4,  1972  Kola  Peninsula 
explosion  as  the  Lg  waves  are  propagated  across  the  continental  path  of  the  Baltic 
shield  only. 

c.  The  wavetrain  of  Lg  is  completely  absent  for  the  December  4,  1988  Novaya 
Zemlya  event  to  ARCESS  and  has  extremely  low  amplitude  for  the  October  25,  1984 
Novaya  Zemlya  event. 

4.  Comparison  of  the  finite-element  simulated  results  of  the  October  25,  1984 
Novaya  Zemlya  explosion  with  the  NORSAR  trace  NOlAO(b)  shows  that  the  time 
of  the  arrival  of  Sn  at  NORSAR  agrees  with  the  finite- modeling  calculation  perfectly. 
The  arrival  of  early  Lg  with  small  amplitude,  which  was  interpreted  by  Baumgardt  as 
a  converted  Lg  at  the  continental  margin  of  the  Baltic  shield,  also  agrees  the  finite- 
element  simulated  trace.  Although  Lg  to  Sn  conversion  and  Lg  scattering  might 
also  have  taken  place  at  the  Baltic  shield  continental  margin,  such  conversion  and 
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scattering  did  not  cause  any  major  changes  in  the  already  small  amplitude  Lg  waves, 
which  was  revealed  by  finite-element  modeling  and  observations. 

On  the  basis  of  finite-element  modeling  and  the  above  three  direct  comparisons 
between  the  finite-element  modeling  results  and  the  observational  data,  the  blockage 
of  Lg  for  the  propagation  paths  from  Novaya  Zemlya  to  ARCESS  and  to  NORSAR 
is  caused  by  the  presence  of  the  low- velocity  sediments  in  the  Barents  Sea. 
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